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First principles molecular dynamics study of amorphous Sii_ x Ge x :H alloys 

T. A. AbtewQ and D. A. Drabokfl 

Department of Physics and Astronomy, Ohio University, Athens, OH 4 5701 
(Dated: February 6, 2008) 

We study the structural, dynamical and electronic properties of amorphous Sii-xGe^H alloys 
using first principles local basis molecular dynamics simulation. The network topology and defects 
of the amorphous network have been analyzed. Structural changes and an increase in number of 
defects have been found as the Ge atomic percentage increases from £=0.1 to x=0.5. The electronic 
properties show a narrower band-gap and increased mid-gap and band-tail defect states in the gap 
of the spectrum as Ge composition increases. Investigation on the band tails of the density of 
states show an exponential (Urbach) behavior. The mobility gap is predicted as a function of Ge 
concentration. 

PACS numbers: 61.43.Bn, 61.43. Dq, 71.20.Bc, 71.23-k 
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I. INTRODUCTION 

Hydrogenated amorphous Si-Ge alloys are important 
materials for uncooled microbolometer applications for 
night "thermal" vision and IR sensingi The materials 
are also of basic interest as a mild form of alloy disorder 
(here "mild" refers to the chemical similarity of the two 
species) juxtaposed with topological disorder. Such al- 
loys have been discussed for photovoltatic applications-, 
and features of the electronic structure in photoemission 
measurements 3 . While empirical models of amorphous 
phases have appeared^, no first principles simulations 
have been reported. Our aim is to offer small but reliable 
models of these important materials, compare them to 
experiment, and to make specific predictions of the struc- 
tural origins of defect states appearing near the Fermi 
level. 

The paper is organized as follows. In Sec.lTTl we discuss 
the approximations and parameters used in the ah initio 
code employed here and describe procedures for generat- 
ing a-Sii-^GeziH alloy models. In Sec. IIIII we describe 
the atomic structural properties by studying partial pair 
correlations, coordination and bond angle distribution of 
the network. The electronic properties of localized mid- 
gap and band-tail states are presented in Sec. IIVI The 
dynamical properties and vibrational density of states are 
given in Sec.0 Finally, we present conclusions. 



II. METHODOLOGY 

A. Total energy and force code 

The density functional calculations in the present 
work were performed within the generalized gradient 



approximation 6 (GGA) using the first principles code 
SIEST A 7 i 8 i 9 . Calculations in this paper employed single 
£ basis with polarization orbitals (SZP) for Si and Ge and 
a double £ polarized basis (DZP) for hydrogeniS. The 
structures were relaxed using conjugate gradient (CG) 
coordinate optimization until the forces on each atom 
were less than 0.02 eV/A. We used a plane wave cut- 
off of 100 Ry for the grid (used for computing multi- 
center matrix elements) with 10 for the tolerance of 
the density matrix in self consistency steps. We solved 
the self-consistent Kohn-Sham equations by direct diago- 
nalization of the Hamiltonian and a conventional mixing 
scheme. The T(k = 0) point was used to sample the 
Brillouin zone in all calculations. 



B. Model Formation 



Realistic models of a-Si have been obtained from the 
WWW algorithm 11 . To model hydrogenated structures, 
we developed a 223 atom a-Si:H model by removing two 
Si atoms and adding nine H atoms (of the 216-atom a-Si) 
to terminate all the dangling bonds except one (to enable 
the observation of one dangling bond defect). Each of the 
H atoms were placed about 1.5 A from the corresponding 
3-fold Si atom. After relaxation, we replaced some of the 
Si atoms by Ge atoms at random, and then relaxed the 
alloy using conjugate gradient minimization to generate 
a-Sii-^GeajiH alloys, with the Ge fraction x being 0.1, 
0.2, 0.3, 0.4 and 0.5. Note that such clustered H agrees 
with the proton NMR second moment data in a-Si:H 12 . 



III. ATOMIC STRUCTURAL PROPERTIES 

A. Pair correlation functions 

The topology of models may be analyzed by partial 
pair correlation functions g a p(i~) of atomic species a and 
(3. The partial pair correlation g a p(r) can be written as 

9*p( r ) = 7— o-T7- : y]S(r-nj) (1) 



4irr 2 pNc a cp 



Mi 



where AT is the total number of particles in the system; 
p = y is the number density , c a = ^ and cp = 
We have used a Gaussian approximation for the delta 
function with broadening a = 0.1 A. 
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FIG. 1: (Color Online) Partial pair distribution function g a /3 
of a-Sii-zGe^iH alloys for two (x—0.1 and x=0.5) Ge atomic 
percentage of (a) Si-Si, (b) Ge—Ge, and (c) Si-Ge. 



We have analyzed the five alloy compositions, and a 
clear first nearest neighbor peak with subsequent deep 



minimum is evident. In Fig.^we plot the partial pair cor- 
relation for Si-Si, Ge-Ge, and Si-Ge in the a-Sii-^Ge;,;:!! 
alloy for £=0.1 and a;=0.5. As the Ge content increases 
from x=0.1 to x=0.5, the peak is fixed at r ^2.50 A. 
The partial correlation for Si-Ge has a peak at r a ^2.4 
A. For the same Ge atomic percentage, the partial pair 
correlation of H-H, Si-H, and Ge-H show a first peak 
r of ^1.45 A, ~1.5 A and ~1.6 A respectively. These 
results reproduce a trend observed in experiment^. 



B. Geometry of Bonding 

We obtained partial ncc ^Sb "h, and average coor- 
dination numbers, n, based on the nearest neighbor dis- 
tance determined in the preceding section, first neighbor 
coordination numbers n S iSi, nsiGc "-siH, n Go Gc, "-GoSi, 
and nceH are obtained by integrating the pair correlation 
function Airr 2 pg a p(r). The results are shown in Table [I] 
The total coordination numbers for Ge, Si are, n Go — 
«GcGc + "GcSi + «-GcH and n Si = n SiSi + n SiGo + n S iH 
respectively. We have obtained 1.0 for total coordination 
of H in all cases. 



TABLE I: The value of r in the first peak of the g(r) and 
the first neighbor coordination number n a /3 in a-Sii-^GexiH 
alloys for five different Ge atomic compositions £=0.1-0.5. 
The integration ranges are from 0.0-2.8 A for Si-Si, Ge-Ge, 
Si-Ge, Ge-Si, and 0.0-1.8 A for Si-H, Ge-H, H-Si, and H-Ge. 
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Si-Si 
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Si-Ge 
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0.06 


Ge-Ge 
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Ge-H 


1.60 


0.04 


0.02 
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H-Si 


1.53 


0.89 


0.89 


0.78 


0.78 


0.67 


H-Ge 


1.60 


0.11 


0.11 


0.22 


0.22 


0.33 



To investigate the effect of Ge composition on the 
structures, we analyzed and obtained all types of bonding 
and structures in the network for each Ge compositions 
considered. For x=0.1, about 96.81% of Si and 95.83% 
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TABLE II: Average percentage m a (T) (bold characters) of 
atoms of species Si and Ge, I fold coordinated at a distance of 
2.68 A for both Si and Ge, and 1.5 A for H in a-Sii-^Ge^H 
alloy for Ge atomic composition £=0.1. We also give the 
identity and the number of Ge and Si neighbors for each value 
of m a (l)- 



TABLE III: Average percentage m a {l) (bold characters) of 
atoms of species Si and Ge, I fold coordinated at a distance of 
2.65 A for both Si and Ge, and 1.5 A for H in a-Sii-aiGe^:!! 
alloy for Ge atomic composition a;=0.5. We also give the 
identity and the number of Ge and Si neighbors for each value 
of m a (l)- 



Si I = 3 1.06 Si 1 = 3 2.32 

Si 3 0.53 Si 2 Ge 1.16 

Si 2 Ge 0.53 Si 3 1.16 
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Ge 1 = 3 4.17 1=4 95.83 

Si 3 4.17 Si 4 79.17 



Si 3 Ge 12.50 Ge 1 = 3 0.96 

Si 2 GeH 4.17 Si 2 Ge 0.96 



of Ge are fourfold, only 1.06% of Si and 4.17% of Ge 
are threefold coordinated, 2.13% of Si are fivefold. No 
fivefold coordination is obtained for Ge. Where fourfold 
Si is concerned, Si4 is a dominant structure which is fol- 
lowed by Si 3 Ge. We observed a similar pattern in the Ge 
fourfold coordination that Ge bonded to S14 structure is 
highly dominant which is followed by Ge bonded with 
Si 3 Ge. The detailed results are shown in Tabic ITT1 

In the case of £=0.5, we observed ~10.47% fivefold 
bonds for Si. About 87.21% of Si and 99.04% of Ge are 
fourfold, only 2.33% of Si and 0.96% of Ge are three- 
fold coordinated. The dominant structure in this case 
is Si atom bonded to Si 2 Ge 2 which is followed by Si 
bonded with Si 3 Ge. Similar pattern is observed in the 
Ge fourfold coordination that Ge bonded to Si 2 Ge 2 is 
the primary structure which is followed by Ge bonded 
with Si 3 Ge. The detailed results are shown in Table ITTT1 
Comparing the bonding statistics with the initial config- 
uration, where we have randomly substituted Si atoms 
by Ge, upon relaxation we have observed: a decrease in 
the number of four folded atoms, an increase in the num- 
ber of strained bond (2.5 A< r <2.7 A), ~ 9.5% for the 
case of x=0.5, and decrease in the normal bonds (r < 
2.5 A). 
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1. Angular Distribution 

We have calculated the partial angular distribution for 
a-Sii_ x Gej;:H with £=0.1, £=0.3, and £=0.5 Ge com- 
positions and plotted in Fig.^Ja)-(f). Though we chose 
only the three Ge compositions the results reported here 
are similar to the other two Ge compositions cc=0.2 and 
£=0.4. The partial pair correlation functions for a-Si-7 
are plotted in the upper panel and the partial pair cor- 
relation functions for a-Ge-7 are plotted in the lower 
panel. In each of the cases considered, we found total 
angular distribution peaks centered near the tetrahedral 
angle with 9 in the range 103°-110°. 

Analysis of different Ge concentrations showed similar 
total angular distribution (the sum of all the partial an- 
gular distributions) because of the similarity in chemical 
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properties and the proximity in different bond-lengths of 
Si-Si, Si-Ge, and Ge-Ge. However, as mentioned above, 
the different compositions affects the type of bonds and 
structures formed in the network. There are also other 
partial angular distributions coming from H. The mean 
partial angular distribution of H-Si-Si and H-Ge-Ge are 
close to the tetrahedral angle 109.47° while the other two 
partials, H-Ge-Si and H-Si-Ge, are off from a tetrahe- 
dral angle and ranges from 100. 0°— 1 16.0°. The probabil- 
ity density for cos (9) is normal (Gaussian). Normally dis- 
tributed cosines of bond angles has been associated with 
exponential band tails in the electron density of states 
near the valence and conduction edg 



IV. ELECTRONIC STRUCTURE 

A. Density of states 

Electronic structure has been described by the elec- 
tronic density of states (EDOS), which was obtained by 
summing suitably broadened Gaussian centered at each 
Kohn-Sham eigenvalue. In the five different Ge composi- 
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FIG. 2: (Color Online) The partial bond-angle distribution 
function as a function of bond angle 8 in a-Sii-xGe^:!! for 
x—0.1 (left panel), for x=0.3 (middle panel), and for x=0.5 
(right panel), (a), (b) and (c) are partial angular distribution 
for three possible angles centering Si and (d), (e), and (f) are 
partial angles taking Ge as a center. 



tions we analyzed the EDOS. The results showed a nar- 
row band-gap spectrum ~1.6 eV for a;=0.1, that becomes 
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FIG. 3: (Color Online) The electron density of states for a- 
Sii-xGea,:!! for £=0.4. The exponential fit in the inset for 
the valence band tail is 2.5 y.e~ E/E ° with E a = 192 meV for 
a;=0.4. The Fermi level is shifted to E=0. 



In this section, we present the results of one of the 
alloys (for a-Sii^Gea,:!! with a;=0.4) among the five dif- 
ferent Ge atomic compositions. The electron density 
of states (EDOS) shown in Fig. shows a narrow gap. 
The band tails of the spectrum which we take in the re- 
gion (-0.7-0.0 eV ) fits exponential with ~ e^ E l E ° with 
E Q = 192 meV as shown in the inset of Fig. 01 We an- 
alyze these defect states in the spectrum with detail in 
the next sections. 



B. Localization: Inverse participation ratio 

In order to understand the electron localization we 
used the inverse participation ratio, IPR, 
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FIG. 4: (Color Online) Inverse participation ratio (IPR) in the 
a-Sii-ajGe^iH alloy for x—0.1 versus the energy eigenvalues. 
The dashed line is the Fermi level. 



where qi(E) is the Mulliken charge residing at an atomic 
site i for an eigenstate with eigenvalue E that satisfies 
[qi(E)] — 1 and N is the total number of atoms in the 
cell. For an ideally localized state, only one atomic site 
contributes all the charge and so IPR=1. For a uniformly 
extended state, the Mulliken charge contribution per site 
is uniform and equals 1/N and so IPR=l/iV. Thus, large 
IPR corresponds to localized states. 

To investigate how the electronic properties change in 
the gap, we have calculated the IPR of a-Sii-^Ge^iH al- 
loy for two different Ge compositions, £=0.1 and £=0.4. 
In each case, we calculated the IPR of the alloy in its 
relaxed ground state. We have also obtained contribu- 
tions of each of the atoms in the alloy to the total IPR 
of a particular eigenstate. This allow us to see different 
types and geometries of structures which are responsi- 
ble for generating a given state. Since we are interested 
in band-tail and mid-gap states, we narrowed our pre- 
sentation here only to those eigenstates which are in the 
mid-gap or near the band-tails of the spectrum. 

We plotted the inverse participation ratio and the con- 
tributions of each of the atoms to the IPR for x=0.1 are 
plotted in Fig.^land Fig. [S] respectively. For the IPR, we 
only chose those atoms which contribute 10% and more 
for a particular state labeled (a-f). Those structures in 
the alloy which correspond to the selected mid-gap and 
band-tail states labeled (a-f) are shown in Fig. EJ As 
we can see from Fig. there is one mid-gap state and 
about five band-tail states. We estimated the mobility 




FIG. 5: (Color Online) The contribution of atoms to the IPR 
(10% and above) of a given state in a-Sii-zCe^:!! alloy for 
£=0.1. The labels from a-f corresponds to different mid-gap 
and band-tail states of Fig. 0] 



band-gap in this case to be (~1.6 eV). 

The structures which are responsible for the mid-gap 
state labeled c, come from a threefold Si, a fivefold Ge, 
and geometrical defects. We have also obtained geometri- 
cal defects corresponding to three of the band-tail states 
labeled (a, d, and e). The larger IPR comes from one 
of these geometrical defects which contains 3 fourfold 
Si atoms connected to each others which have angular 
distribution off from a tetrahedral angle by ±15°. The 
larger contributions to the other two mid-gap states come 
from a fivefold Si atom together with a geometrical de- 
fect (6), and a threefold Si bonded with a fivefold Si and 
a geometrical defect (/). 

In the case of x=0A, the inverse participation ratio as a 
function of eigenvalue is plotted in Fig. [7| while the con- 
tributions of each of the atoms to the IPR (only those 
atoms which contribute 5% and more) for a particular 
state labeled {a-f) and the different structures associ- 
ated with these states are plotted in Fig. |S| and Fig. 
respectively. As we can see from the IPR plots in the two 
cases, ir=0.1 and x=QA, as the Ge content increases, we 
observed an increasing number of band-tail states close to 
the conduction band edge and hence a narrow bandgap 
spectrum. 

Further analysis of the type of bonds and structures 
associated with the mid-gap states and band-tail states 
showed that these states are coming not only from the 




FIG. 6: (Color Online) Structures which contributes 10% and 
above to the inverse participation ratio that corresponds to 
the mid-gap and band-tail states shown in Fig.|l|from a-f in 
the a-Sii-^Ge^iH alloy for x=0.1. The color codes are blue 
for Si atoms and red for Ge atoms. 



dangling bonds defects (threefold and twofold bonds) but 
also from geometrical defects where there exist a strain 
bond angle in the structure which has an angular dis- 
tribution far (±15°) from a tetrahedral angle. As com- 
pared with the case where x=0.1, we have observed an 
increase in the contribution from the geometrical defects 
to band-tail states for x=0A. As shown in Fig. [§1 the 
larger contributions to the state labeled (a and b) come 
from geometrical defects, a fivefold Si, and a threefold 
Ge structures. States labeled d and / in this case are 
dominated by a geometrical defect which has 3 four- 
fold Si atoms connected to each other with strain bond 
with angular distribution off from a tetrahedral angle by 
±20°. The number of Ge atoms in these geometrically 
defected structures has increased as compared with the 
case £=0.1. The dominant contributions to the mid-gap 
state arise from a fivefold Si atom bonded with two geo- 



0.06 



Q_ 




Eigenvalues (eV) 



FIG. 7: (Color Online) Inverse participation ratio (IPR) 
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FIG. 8: (Color Online) The contribution of atoms to the 
IPR (5% and above) of a given state in a-Sii-xGe^H alloy 
for :r=0.4. The labels from a-f corresponds to different mid- 
gap and band-tail states of Fig. 



metrical defects. 

To emphasize the effect of Ge concentration in the mo- 
bility band gap, we have plotted the mobility gap as a 
function of the Ge concentration x in Fig.E3(f)- The mo- 
bility gaps are extracted from Fig.^|(a)-(e) which show 
the inverse participation ratio as a function of eigenval- 
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FIG. 9: (Color Online) Structures corresponding to the mid- 
gap and band-tail states shown in Fig. |7| from a-f in the 
a-Sii-zCe^R alloy for x—0.4. The color codes are Si atoms 
(blue) and Ge atoms (red). 



ues for different x. We observed a decrease in the mobil- 
ity gap as the Ge concentration increase from £=0.1 to 
x=0.5. The mobility gap is estimated by examination of 
the plots of IPR versus energy. In each case, there is a 
reasonable well-defined energy near the valence and con- 
duction tails at which the IPR becomes roughly constant 
reflecting extended states. We include "error bars" to 
convey a rough estimate of uncertainty in our estimated 
gaps. 
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FIG. 10: Inverse participation ratio (IPR) versus the energy 
eigenvalues for different x: (a)-(e) and the estimated mobil- 
ity gap plotted versus the Ge concentration x: (f), in the 
a-Sii-xGe^H alloy. 



ing each atom with 0.03 A in three orthogonal directions 
and then performing first principles force calculations for 
all the atoms for such displacement. Each calculation 
yields a column of force constant matrix. The vibra- 
tional eigenvectors and eigenvalues of the supercell are 
then obtained by diagonalizing the dynamical matrix. 

In Fig. ^2 the phonon density of state for a- 
Sii-^Ge^H for a;=0.4 is plotted. The acoustic peak ap- 
pears at 13.29 meV and two optical peaks appear at 31.98 
meV and 49.01 meV which is in good agreement with the 
experimental result reported by Mackenzie et a<i£. The 
higher frequency modes in the range (213 meV-236 meV) 
are associated with hydrogen atoms with H-Si and H-Ge 
bonds which is in agreement with the experimental result 
reported by Wells et ali&. The principal hydrogen re- 
lated features of the spectrum which exhibit higher IPR 
(highly localized states) are: stretch modes of Si-H at 
2039 cm -1 , 2013 cm" 1 , 1980 cm" 1 , and 1873 cm" 1 , and 
of Ge-H at 1904 cm -1 , 1646 cm -1 and a wagging vibra- 
tion modes of Si-H and Ge-H dominates in the region of 
600 - 900 cm" 1 . 



V. DYNAMICAL PROPERTIES 

The dynamical properties of a-Sii-zGe^H alloys are 
analyzed with the vibrational density of states (VDOS) 
and inverse participation ratio. The vibrational eigen- 
values and eigenvectors are obtained from the dynamical 
matrix. The dynamical matrix is determined by displac- 



CONCLUSION 

We have presented an ab initio first principles molecu- 
lar dynamics study of a-Sii-^Ge^H alloys for five differ- 
ent Ge atomic compositions. In this study, we reported 
a systematic investigations of structural and electronic 
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FIG. 11: (Color Online) Vibrational density of states (black) 
and the inverse participation ratio (blue) for the a-Sii-^Ge^iH 
model for £=0.4. 



properties of a-Sii-^Ge^iH alloys and the consequences 
of Ge atomic composition in the structure and electronic 
properties of the alloy. Where the structural properties of 
a-Sii-^GeiciH alloys are concerned our results definitively 
show (a) very small change in the total coordination and 
total bond angle distribution, and (b) the emergence of 
preferential structures in the network as the Ge atomic 
content increases. The electronic density of states shows 
an increase in the band-tail states and narrow band gap 
as the Ge content increases due to the increase in the 
number of defects. This suggest the potential of Ge (by 
varying its atomic compositions) in tuning the band gap 
of the &-Sii- x Ge x :ll alloy for different technological ap- 
plications. 
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